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Abstract. The availability of Bose-Einstein condensates as mesoscopic or macro- 
scopic quantum objects has aroused new interest in the possiblity of making and de- 
tecting coherent superpositions involving many atoms. In this article we show that it 
may be possible to generate such a superposition state in a reasonably short time using 
Feshbach resonances to tune the inter-atomic interactions in a double-well condensate. 

We also consider the important problem of distinguishing whether a coherent 
superposition or a statistical mixture is generated by a given experimental procedure. 
We find that unambiguously distinguishing even a perfect 'NOON' state from 
a statistical mixture using standard detection methods will present experimental 
difficulties. 



1. Introduction 

The idea of a macroscopic state demonstrating quantum mechanical behaviour was 
introduced by Schrodinger in 1935 [I]. His famous thought experiment considered how 
a macroscopic entity (in this case a domestic cat) could evolve into a superposition 
of two distinct physical states (alive and dead) when entangled with a microscopic 
system that obeyed the laws of quantum mechanics. This possibility of a macroscopic 
system being simultaneously in two distinct physical states was initially considered to 
be a flaw in quantum mechanical theory pQ. However, experiments have shown the 
predictions of quantum mechanics to be correct as superpositions of macroscopically 
distinct physical states have been produced in a variety of systems. These are often 
referred to as 'Schrodinger cat states'. For example, a group of six beryllium ions has 
been put into a superposition of two hyperfine states [2], and persistent currents (of a 
few fiA) of opposing circulation in SQUIDs have been detected [3]. 

Bose-Einstein condensates (BECs) are attractive candidates for generating 
macroscopic superpositions due to the large number of atoms that share a single 
quantum state. This may provide a useful system in which to further test the 
validity or boundaries of the assumption of macroscopic realism. A macroscopic 
superposition has yet to be demonstrated in a Bose-Einstein condensate, although 
there exist numerous proposals for generating either a superposition of relative phase 
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or number states [SI El El El El EED] • In the following we consider a particular kind of 
superposition state of a single component BEC in a double well, that is, a superposition 
of the two states where the entire condensate is localised in one of the wells (sometimes 
referred to as a 'NOON' state). If realisable, this kind of superposition promises to 
be useful in quantum information applications and precision interferometry, due to the 
measurement uncertainty scaling inversely with the number of particles (the so-called 
Heisenberg limit). 

A major difficulty in realising macroscopic quantum superpositions is decoherence, 
which occurs when interactions with the surrounding environment cause the pure 
superposition state to decay into a statistical mixture. However, this paper 
is not concerned with avoiding decoherence in the realisation of macroscopic 
superpositions. Instead, we concentrate on measurements aimed at distinguishing a 
coherent superposition from a statistical mixture. We find even in a decoherence-free 
environment, demonstrating a superposition presents several practical challenges. Once 
realised, measurements of the purity of the state could be useful in studying, for example, 
rates of decoherence. 

This paper begins with an introduction of the two-mode description of a double- well 
condensate in Section [21 In Section [3] we consider measurements aimed at distinguishing 
between a coherent superposition and a statistical mixture. We focus on quadrature- 
based measurements, analogous to those used in quantum optics, that can be realized 
with a Ramsey-type interference experiment. The coherence of the NOON state is 
evident in parity measurements of the number distribution after the two modes are 
interfered — a difficult measurement with standard atom counting techniques. We 
analyse the effects of atomic interactions and discuss atom loss during the interference 
procedure. Finally, we show in Section H] that a mesoscopic superposition of 20 atoms 
could be generated in a reasonable time frame if a Feshbach resonance can be used to 
tune the atomic interactions, before concluding in Section [51 



2. Theoretical Models 



2.1. Two-mode approximation 

The Hamiltonian for a condensate in an external trapping potential, V ex t(r), is 
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where ip is the field operator for the condensate, and the non-linear interaction parameter 
is Uq = Awah/ m, (a is the s-wave scattering length describing two-body collisions within 
the condensate, m is the atomic mass). We consider the case where the external potential 
provides a double well confinement for the condensate. Double well potentials can 
be generated by an optical lattice with an additional harmonic confinement to reduce 
the number of occupied lattice sites to two pTj . They can also be realised on chips, 
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where suitably arranged current carrying wires create a magnetic confinement for the 
condensate atoms |12j . 

When a double well potential is considered, the above Hamiltonian can be simplified 
by making use of the two-mode approximation. This means we consider each atom to 
be in some linear superposition of being in the left well and being in the right well. 

We consider the zero temperature case, where all atoms in the system are condensed. 
If the ground state energies of the condensate in the two single (and separate) wells are 
sufficiently separated from the energies of the condensate in all other excited single 
particle states, transitions to or from the two modes of interest and these higher lying 
states can be neglected. This is required for the two-mode description to be valid. In 
the two-mode approximation, the field operator is expanded as 

ip{r) » (p L (r)a L + <p R (r)a R , (2) 

where and a R are discrete Bose annihilation operators for the left and right well 
respectively, and 4>l/r are the ground state spatial wave functions of the condensate in 
the left and right wells. 

Substituting this into equation <^f, we find an effective Hamiltonian 

H ef f = hE L a[a L + HE R a R a R + -^-a[a[a L d L (3) 

+ —^-a R a R a R a R - Kk [a L a R + a R a L j , 

where we have neglected the spatial overlap of the left and right well densities. The 
single well bound state energies, El/ r , are 

El/r = \J " dr p m lf) (^V 2 + V ext (r)^j 0^(r). (4) 

k, the tunnel coupling, is 

k = df<Pl /R (r) (^V 2 + VUA (5) 

and the effective non-linear interaction terms are 

U L/R = U J dr\<p L/R (r)\\ (6) 

For the remainder of this paper we assume a symmetric potential, where El = E R = 
and U L = U R = U. 



2.2. Fixed number representation 

By ignoring all possibility of atom loss and decoherence, we can efficiently represent the 
iV-body wave function using the basis \N — n,n), representing states with N — n atoms 
in the left well, and n atoms in the right well. Any wave function can be written as a 
superposition of these number states, i.e. 

N 

m)) = Y. c ^\ N - n ' n ^ ( ? ) 

71=0 
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where Y2 n \ c n(t)\ 2 = 1- in this representation, the expectation value of the number of 
atoms in the left well is 

(N L ) = J2( N - n )\cn(t)\\ (8) 

n 

and the variance in the number difference is 



V = ^{N -2nf\c ri 



J2(N-2n)\c r , 



(9) 



For any initial state, these coefficients have a time-dependence given by 

c n (t) = {N-n, n\ip(t)) = (N-n, n|e-^ /ft |^(0)). (10) 
Taking the derivative with respect to time, we find [13] 



^M = -l(N-n,n\n\m)- (11) 

Inserting the Hamiltonian, (JTJ, into this expression gives the equations of motion for 
the number state coefficients 

ih 6 ^- = hU L [(N - nf -{N- n)] c n (t) + KU R [n 2 - n] c n (t) 

- Hk^(N -n + l)n c n+1 (t) - hn^J(n + 1) (N - n) c n _ 1 (t)(12) 

The ground states and dynamics of a condensate in the double well can then be found 
by solving these equations numerically. 

Fig. [1] shows the ground state probability distributions for the atom number in 
the left well (found using imaginary time propagation p3]) for three different regimes 
- negligible, intermediate, and strongly attractive interactions. The variance in the 
number difference for these are 12, 293, and 396 respectively (a perfect NOON state 
would have a variance of 400). 

When non-linear interactions are weak compared to the tunneling, the ground 
state number distribution is essentially binomial (Fig. [T] (a)). For large N this can 
be approximated as a Poissonian distribution and the ground state would be a coherent 
state (which has equal mean and variance of atom number in each well). For intermediate 
attractive interactions, ground states exist consisting of two well separated peaks in the 
number state distribution. These can be considered as superpositions of distinct physical 
states, and we refer to them as 'mesoscopic' superpositions 

The ground state in the limit of infinitely strong attractive interactions is a 
macroscopic superposition of the entire condensate localized in each well. This can 
be understood by realizing that it is energetically preferable for the condensate to be 
localized in one well, but that given the symmetry of the double well potential these 
two localized states are degenerate. Thus in the ground state it is equally likely that 
the condensate will be found entirely in one well as in the other (hence the two peaks 
in the probability distribution). This state is the macroscopic superposition that we are 
interested in. Fig. [1] (c) shows that with a ratio of non-linear interaction strength to 
tunnelling rate U/k— —0.5, the ground state of the double well condensate is close to, 
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Figure 1. Ground state probability distributions for weak, intermediate, and strong 
attractive interactions. Nl is the atom number in the left well. Interaction strengths 
are (a) U/n = 0, (b) -0.1, and (c) -0.5. 



but not quite, an ideal superposition. We return to the problem of creating such a state 
in Sec. HJ 

3. Verifying a superposition state 

We now consider the problem of demonstrating that a state is a macroscopic quantum 
superposition, as opposed to a statistical mixture. This will allow tests of macroscopic 
realism, may be an important tool for studying decoherence, and is a prerequisite to 
building a practical interferometer using NOON states. 

The obvious measurement to make on such a system is the atom number in each 
well. For the NOON state, we have equal probability of finding N atoms in the left well 
and none in the right, or vice- versa. However, such a measurement can not distinguish 
the coherent NOON state 

|0) = (|iV,O) + e^|O,iV))/v / 2 (13) 

from the statistic mixture with density operator 

/3=(|iV,0)(iV,0| + |0,iV)(0,iV|)/2. (14) 

There exist complimentary measurements that distinguish between a coherent 
superposition and a statistical mixture. One method of finding such an measurement 
is motivated by the expansion of the density operator of the statistical mixture as the 
average of all coherent phases 4>, 

We see that the determination of constitutes proof of coherence — as is undefined 
for the mixed state. 

Consider an arbitrary operator x; for the mixed state in (fl4|) . the expectation value 
of x is the average of the two pure, separable states, 

(x)mixed = \(( N , 0) + (0> N\x\0, AO) • (16) 
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On the other hand, for the pure NOON state the expectation value is 




(17) 



The additional interference terms clearly display a dependence on the value of 0. 
Therefore, an appropriate observable x has non-zero (0, N\x\N, 0). An example of an 
operator that achieves this is d^d^ + d^d 1 ^, which coherently transfers N atoms from 
one well to the other or vice-versa. Unfortunately, there is no clear way of directly 
measuring this observable in an experimental setting. 

3.1. Quadrature phase measurements 

We now consider a quadrature-based method for distinguishing the entangled NOON 
state from the statistical mixture. In this paper, we define the quadrature operator as 



The measurement of this observable can be achieved using simple linear interference 
and number measurements. Such a procedure is analogous to quantum optics 
experiments using a 50-50 beam splitter to interfere two photonic modes before intensity 
measurement, and allows access to phase information. Unlike common quantum optics 
experiments, both modes contain a similar number of atoms and neither mode can be 
interpreted as a local oscillator. This explains the difference between the above definition 
and the standard quadrature arising from homodyne measurement, proportional to 
die 19 + d\e~ %e . Both definitions have been employed in theoretical discussions of BECs 
in the past pi [HUE]. 

To realize these quadrature measurements, we propose a Ramsey-type experiment 
(see Ref. [20] for a description of the Ramsey technique). Similar experiments have been 
proposed for double well condensates (for example, to detect a weak force [21]). The 
first step, after creating the superposition state, is to set the tunneling rate between 
the wells n and interaction strength U to zero and let the system evolve for some time, 
5t, during which an energy imbalance exists between the wells (i.e. El — E^ = 5E is 
non-zero). The quadrature angle is set by 6E5t = 9. The second stage is to restore the 
symmetry of the wells (i.e. set 5E to zero), and switch on tunneling for a time of 7t/4k 
(this is analogous to a beam splitting operation). After this we make a measurement of 
the atom number in each well. The difference in atom number is exactly proportional to 
Xg. We simulate the entire procedure using the equations of motion of the number state 
coefficients ( 1121) . It is then straightforward to numerically calculate the distribution of 
measurement outcomes, as given by |c n | 2 , and thus any moment (e.g. mean, variance, 
etc) of the the quadrature Xg. 

The next step is to extract information about the off-diagonal terms of the density 
matrix by measuring the interference terms in ( fl7l) . It is straightforward to see that 
(N,0\Xg\0, N) = for N > 2. In the special case of N = 1, interference is observed 
by a sinusoidal dependence on the value of 9 — <p, as seen in Fig. [2] (a) (where is the 



X g = d L d R e 10 + d R d L e %e . 



(18) 
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Figure 2. (a) Interference fringes of (Xg) as a function of the accumulated phase 
shift for an ideal superposition containing 1 atom, (b) The same quantity for an ideal 
superposition containing 2 atoms displays no interference fringes, (c) For the two atom 
superposition state, the variance (V(Xg) = (Xg) — (Xg) 2 ) is sensitive to the phase 9 
accumulated during the Ramsey simulation. 



phase angle used in (13)). In Fig. [2] (b) we plot the same quantity for N = 2, and it is 
unsurprising that inteference is lacking for this case. 

However, higher moments of the quadrature measurements do contain information 
that can distinguish a pure NOON state from a classical mixture. For N atoms, the iVth 
moment (Xg*) contains exactly one subterm equal to a L + a^a^ , while all lower 
moments lack such a term. In Fig. [2] (c), we see interference fringes in the quadrature 
variance for the case N = 2. The frequency of these fringes is doubled compared to the 
case N — 1. 

In general the N atom NOON state will display fringes with frequency N times 
greater than the first-order coherence measured by (Xg) allows. Specifically, the iVth 
quadrature moment contains terms proportional to cos(iV(# — 0)). It is this scaling 
that makes NOON states of interest for precision interferometry — a NOON state with 
known <ft could potentially be used to measure SESt with accuracy proportional to N~ l 
(the so-called Heisenberg limit), compared with the A r_1 / 2 scaling typical when using 
'classical' inteferemetric techniques. Such scaling has been observed in single-photon 
experiments [23], but to-date neither with atoms nor NOON states. 

The quadrature moments are intrinsically linked with the number distribution \c n \ 2 
after the Ramsey interference procedure, and so it follows that the off-diagonal terms 
in Eq. (1X7]) are directly visible in this distribution. We have plotted the output in 
Fig. [3] (a) for 9 = n/2, and we observe a pattern where each second |c n | 2 is zero. This 
interference pattern is sensitive to the accumulated relative phase, so that for certain 
values of 9 — <p the interference pattern is absent and the number state coefficients are 
given by a binomial distribution. On the other hand, if the initial state were a statistical 
mixture a binomial distribution would be expected for all values of the relative phase. 

The interference can be seen most clearly in the parity, defined as 

p = ^(-iri Cn | 2 , ( 19 ) 
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Figure 3. Atom number distribution |c„| 2 after the Ramsey procedure is simulated, 
with accumulated phase ir/2. Initial state has 20 atoms and <fr — 0. Non-linear strength 
U/k is (a) zero, (b) -0.01, and (c) -0.025. 
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Figure 4. Expectation value of parity after Ramsey simulation using ideal 
superposition state containing (a) 10, and (b) 20 atoms. Non-linear interactions reduce 
the amplitude of the parity oscillations. (Note that in (a) the U/k = —0.01 line lies 
under the U/k = line.) 



which is the difference in probability that an even or odd number of atoms in measured 
in one of the modes at the end of the Ramsey experiment. The expectation value 
of the parity is an oscillatory function of the relative phase between the wells, again 
with a frequency proportional to the total atom number N, as seen in Fig. H] (a,b) . 
Measurements of the parity were considered in detail in Ref. [21], where the authors 
discuss methods for obtaining Heisenberg-limited phase resolution. 

3.2. Measurement difficulties 

Accurate measurement of the iVth quadrature moment, or equivalently the atom parity 
after the Ramsey experiment, is experimentally very difficult. Measuring the interference 
pattern in the c n coefficients would require very accurate atom counting, as the relative 
phase value giving maximum fringe visibility for an even total atom number gives 
minimum fringe visibility if the total atom number is odd, and vice versa. Indeed, it is 
apparent that the loss of a single atom would completely destroy the coherence of the 
NOON state. Therefore, the required counting efficiency rj is such that which 
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would be challenging for large values of N. An alternative method of determining the 
parity has been suggested in Ref . [23] . This method does not require extremely accurate 
atom counting; however it does involve a third condensate mode and precise control of 
the nonlinearity. 

There are additional complications that may make this procedure difficult, even 
for moderate values of N. We now consider the effect of finite non-linear interactions 
during the Ramsey interference procedure. With an ideal superposition as the initial 
state, the non-linear interactions degrade the visibility of the interference pattern in 
the c n coefficients, as shown in Fig. [3] (b,c). The oscillations in the expectation value 
of parity are also reduced by the presence of non-linear interactions (see Fig. H|). For 
a given value of U, the visibility of the parity oscillations decreases as the number of 
atoms in the superposition increase. However, we note that the phase of the interference 
fringes is unaffected by the non-linearity, which is important for possible interferometric 
applications. 

Finally, when the initial state is an ideal NOON state, we observed that including 
non-linear interactions can result in small fringes in the mean quadrature value even for 
N > 2. The maximum amplitude of these induced fringes is much less than the total 
atom number. These fringes are most noticeable for small atom numbers, as they have 
amplitudes of only one or two atoms regardless of the total atom number. 

4. Generating a superposition state 

There have been a number of proposals for generating superpositions of various kinds 
involving Bose-Einstein condensates. Many of these consider superpositions of two- 
component condensates. In this case, methods involving adiabatic manipulation [5], 
and dynamical evolution (making use of the interplay between nonlinear interactions 
and tunneling, with a specific initial relative phase) [6], [71 [8], have been suggested to 
generate superpositions. Superpositions of phase states have also been considered [91 [TO]. 

We are interested in superpositions of number states in a single component 
condensate in a double well potential. One proposal for generating superpositions of this 
kind involves using a Feshbach resonance to produce a sudden change in the interaction 
strength, initiating a dynamical evolution where a macroscopic superposition emerges 
periodically [T5] . 

In Ref. [15], it was claimed that the adiabatic method was not feasible due to the 
infinitely long evolution time required because of the near degeneracy of the ground 
and first excited states in the strongly attractive regime. We find that with realistic 
parameters a superposition of two well-separated wave packets can be generated on a 
time scale of seconds using a smooth change in interaction strength. Although this 
evolution is not necessarily adiabatic it can generate a 'mesoscopic' superposition state. 

To estimate some suitable parameters for our effective Hamiltonian we have 
considered a small condensate of Rubidium atoms. Around a magnetic field strength 
of 155 G there is a Feshbach resonance for 85 Rb atoms, allowing the s-wave scattering 
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Figure 5. Evolution of wave function into a mesoscopic superposition caused by a 
linear decrease in interaction strength from U=l s^ 1 to —3 s" 1 in (a) 0.5, and (b) 4 
seconds. Initial wavefunction is the ground state at U=l s _1 and k = 10 s _1 . 

length to be tuned from at least 2000 a to —200 a (a being the Bohr radius). Assuming 
a cigar-shaped condensate confined by trapping frequencies of 1 kHz in both tight 
directions and 100 Hz in the longitudinal direction, and using a Gaussian approximation 
for the condensate wavefunction, we estimate that the interaction parameter, U, given by 
equation (jHJ) could range from approximately 30 s _1 to —3 s" 1 . We look at simulations 
which are a few seconds in length (note that condensate lifetimes of greater than 10 
seconds are experimentally achievable [16J). 

Fig. [5] shows the wave function as the interaction strength is changed linearly from 
1 s -1 to —3 s _1 over a timespan of 0.5 and 4 seconds. The initial state is the ground state 
at U=l s" 1 , and k = 10 s" 1 . Clearly the slower change in interaction strength results, 
as per the adiabatic theorem, in a final state that is closer to the ideal macroscopic 
superposition. This can be seen in Fig. El which plots the fidelity of the evolving wave 
function, i.e. the overlap, | {(Pnoon^^)) | 2 , of the ideal superposition state with the 
evolving wave function. The variance in the number difference between the wells of the 
final states in Fig. [5] are approximately 283 and 371, compared with 400 for an ideal 
superposition state containing 20 atoms. Once the interaction strength is held constant 
at —3 s -1 the probability distributions \c n \ 2 do not change significantly with time. 

The parameters could be optimized further to obtain better superposition states. 
In general, better states are obtained when the final ratio U / k is large and negative. On 
the other hand, the timescale for the evolution to be adiabatic is inversely proportional 
to both U and n. The requirement of physically separated wells implies a relatively 
small value for k, while the timespan of the experiment is limited by atomic loss. We 
have not attempted to optimize the process of generating a superposition state within 
these bounds, but rather concentrate on measurements aimed at distinguishing even an 
imperfect superposition from a statistical mixture. 

4-1. Intereference measurements of non-ideal superposition states 

We now perform the Ramsey procedure followed by parity measurements in order 
to detect coherence in the non-ideal superposition states generated above. Figs. [7] 
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Figure 6. Fidelity of final wave function compared to the ideal NOON state, for 
different ramp times (dots). Dashed line is fidelity of ground state U = —3 s _1 and 
k = 10 compared to the ideal NOON state. A slower ramp time results in a final wave 
function closer to that of an ideal superposition state. 



and M show typical fringes in the number difference and parity for an initial non- 
ideal superposition state of a 20 atom condensate. In the parity, we see clear high 
frequency components in the interference fringes that correspond to the off-diagonal 
coherence between states of large atom number difference. However, as the state is 
not perfect, other frequency components are present, resulting in beating and a low 
frequency envelope. 

In the Appendix, we formalise the relationship between the parity frequency 
components and coherence and show that 

P = B riim e imrp (N - n, n\p\N -n-m,n + m) + H.c. (20) 

n,m 

for some real numbers B n>m . The component with angular frequency m corresponds to 
coherence between elements separated by m atoms in Fock space. In fact, the highest 
frequency component with angular frequency N is proportional only to (N, 0|/i|0, N), 
and is the only component observed in the results for an ideal superposition (Fig. @|. 
The presence of the same high frequency component in Fig. [8] (a) is an unambiguous 
demonstration that the generated state contains coherence and is not a statistical 
mixture. 

With the introduction of non-linear interactions during the Ramsey interference, 
the amplitude and frequency of the parity fringes become less regular and are no longer 
periodic over 9 — <p modulo 2n. It would not be easy to use the results when U is large 
as an indicator of the presence of a superposition due to their irregularity, and the fact 
that high frequency components cannot be said to indicate NOON-type coherence. 

5. Conclusions 

Our simulations indicate that a quantum superposition of a small Bose-Einstein 
condensate may be generated in a reasonably short time by a smooth change in the 
atomic interaction strength from repulsive to attractive. However, we find that verifying 
that a superposition has been generated may not be straightforward, even in the ideal 
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Figure 7. Expectation value of atom number in each well after Ramsey simulation 
for a 20 atom condensate, for non-linear interactions being (a) zero, (b) —0.1 s^ 1 , and 
(c) —0.25 s~ x . Initial state was generated from the ground state at U=l s _1 , K = 10 
s -1 , and linearly changing the interaction strength to —1 s _1 over 4 seconds. 
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Figure 8. Expectation value of parity after Ramsey simulation of a 20 atom 
condensate, for non-linear interactions being (a) zero, (b) -0.1 s" 1 , and (c) -0.25 s" 1 . 
Initial state was generated from the ground state at U=l s~ x , K — 10 s" 1 , and linearly 
changing the interaction strength to —1 s" 1 over 4 seconds. 



case of a perfect 'NOON' state superposition and no decoherence. We have considered a 
Ramsey-type interference experiment as a method to measure the quadrature operator 
and parity thereof of a double well wave function. 

For a NOON superposition, only the highest-order quadrature moment (and not 
the mean) is sensitive to the accumulated relative phase and is therefore useful in 
distinguishing between a statistical mixture and a coherent superposition. Parity 
measurements after a Ramsey-type experiment on an ideal superposition state show a 
dependence on the relative phase that could in principle be used to verify the presence of 
the superposition. A statistical mixture would not be expected to show any dependence 
on the relative phase if parity measurements were performed. The presence of non- 
linear interactions during the Ramsey experiment degrades the amplitude of the parity 
oscillations. 

Measurements of a non-ideal superposition state also show a high-frequency 
dependence on the relative phase — a smoking gun indication of NOON-type coherence. 
However, we observed that the imperfect state displays additional frequencies resulting 
in a beating pattern. For some choices of phase angle, the visibility is close to one 
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while the period decreases to 27r/N, allowing it to remain useful for Heisenberg-limited 
phase measurements. Similarly as for the ideal state, non-linear interactions degrade the 
amplitude and regularity of the parity oscillations, however for small non-linearity, the 
high frequency component could still be used to verify that a mesoscopic superposition 
has been generated. 

There are many factors that may make these measurements difficult in practice. 
For example, phase diffusion due to non-linear interactions dramatically reduces the 
visibility of the parity oscillations. Other effects not considered here would also reduce 
the fringe visibility, such as 3-body loss, and other sources of decoherence. Highly 
accurate atom counting is required to observe the interference patterns in the c n 
coefficients, and the associated parity oscillations as a function of accumulated phase, 
unless an alternative measurement scheme could be realised. 
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For a state with precisely N atoms, the (discrete) quadrature operator can take N + 1 
values. Therefore, it is possible to express any linear combination of the measurement 
probabilities (i.e. the \c n \ 2 after the Ramsey interference) as a linear sum of the first N 
quadrature moments. Specifically, the parity can be expressed this way, i.e. 

N 

n=0 
N 

= J2M(e i6 a{a R + e- i6 a L aiy), (.1) 

n=0 

where {A n } are real numbers. Expanding the power results in 

N n 

P{0) = E A « m Cne {n - 2m)d (: a[ n - m a^ R m a n R - m : sym ), (.2) 

n=0 m=0 

where : x : sym represents the symmetric ordering of x. 
The operator : 

^ L ^ L ^ R ^ R 

~ m -sym transfers n — 2m = M atoms from the right 
well to the left, and therefore we can say 

(: a[ n - m a™a R m a R ~ m : sym ) =^D n (N- n,n\p\N - n -m,n + m), (.3) 

n 

where D n are appropriately chosen constants. 

Putting everything together, we finally arrive at Eq. (TSUI) , where the the constants 
B n>m can be obtained from the above. 
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